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SUMMARY 


Based  on  the  Rice- Wang  thermodynamic  model  of  intergranular  embrittlement,  total 
energy  calculations  employing  the  IT-APW,  DMol  and  DVM  techniques  address  the  electronic 
basis  of  the  segregation  energy  difference  between  grain  boundaries  (GB)  and  free  surfaces  (FS) 
controlling  the  embrittlement  potency  of  a  segregating  solute.  Precise  FLAPW  calculations 
employing  large  supercells  with  structural  relaxations  determined  by  DMol  predict  the 
embrittling  effect  of  P  in  Fe,  and  indicate  that  B  is  neutral  or  weakly  cohesion  enhancing  in  Fe. 
The  embrittlers  P  and  S  are  found  to  behave  as  "embedded"  atoms  with  negligible  hybridization 
with  Fe,  while  B  hybridizes  to  form  directional  o  bonds  which  significantly  raise  its  relative 
energy  on  the  free  surface.  Magnetic  contributions  are  found  to  decrease  both  the  embrittling 
potency  of  P  and  the  cohesion  enhancing  ability  of  B.  Toward  greater  precision  of  segregation 
energy  calculations,  DMol  calculations  have  been  initiated  with  modified  boundary  conditions  to 
allow  for  GB  free  volume.  DVM  cluster  calculations  are  applied  to  explore  the  third-element 
interaction  with  substitutional  Mo  at  the  GB  core.  A  cohesion  enhancing  effect  of  Mo  more  than 
compensates  the  embrittling  effect  of  P.  GB  fracture  strength  is  measured  in  UHS  steels  and 
correlated  with  segregation  of  P  and  S. 


1.  INTRODUCTION 


The  objectives  of  this  research  are  to  establish  the  electronic-level  mechanisms  of 
intergranular  embrittlement  in  Fe-base  alloys,  and  to  design  and  produce  improved  grain 
boundary  compositions  for  enhanced  stress-corrosion  resistance  in  ultrahigh-strength  steels. 
Based  on  the  Rice- Wang  thermodynamic  model  of  impurity-induced  embrittlement,  our 
electronic  calculations  focus  on  the  origins  of  the  energy  difference  of  solutes  in  the  Fe  grain 
boundary  (GB)  and  free  surface  (FS)  environments  as  the  key  quantity  governing  embrittlement 
potency. 

Our  paper  published  last  year  (Wu,  Freeman  and  Olson  [1])  reported  the  application  of 
the  precise  Full-Potential  Linear  Augmented  Plane-Wave  (FLAPW)  Method  to  comparison  of 
the  energy  and  electronic  structure  of  P  in  the  Fe  L3  [1 10]  (111)  GB  and  corresponding  Fe  (1 1 1) 
FS  environments.  It  was  shown  that  if  Fe  structural  relaxations  are  neglected,  the  chemical 
bonding  of  P  is  actually  greater  in  the  GB  than  the  FS,  in  opposition  to  its  embrittling  effect.  Our 
recently  published  paper  comparing  P  and  B  on  the  Fe  (111)  FS  (Wu,  Freeman,  and  Olson  [2]) 
showed  that  the  degree  of  structural  relaxation  in  the  FS  environment  is  relatively  small.  Also 
recently  published  is  our  study  of  the  detailed  structural  relaxations  of  the  Fe  L3  [110]  (111)  GB 
with  and  without  P,  employing  the  DMol  cluster  method  (Tang,  Freeman,  and  Olson  [3]).  This 
study  showed  that  the  cluster  size  had  to  be  increased  from  53  to  91  atoms  to  fully  treat  the 
relatively  long-range  relaxations  associated  with  the  GB.  The  clean  GB  was  found  to  relax 
locally  to  an  co-phase  structure,  and  a  substantial  energy  penalty  is  associated  with  reconstruction 
of  the  GB  by  P  to  an  "anti-co"  structure. 

The  detailed  reconstruction  of  the  GB  has  been  taken  into  account  in  precise  FLAPW 
calculations  comparing  P  in  the  GB  and  FS  environments  employing  larger  24  layer  supercells 
for  the  GB  and  13  layer  supercells  for  the  FS  in  order  to  achieve  the  most  precise  energy 
comparisons  currently  possible.  A  paper  detailing  these  definitive  calculations  demonstrating  the 
embrittling  behavior  of  P  (Wu,  Freeman  and  Olson  [4])  was  provided  with  our  previous  progress 
report.  A  new  paper  detailing  the  comparison  of  precise  FLAPW  calculations  for  P  and  B  in  the 
GB  and  FS  environments  (Wu,  Freeman  and  Olson  [5])  is  attached  as  Appendix  A.  A  revised 
manuscript  describing  further  DMol  studies  of  the  structural  relaxations  and  bonding  behavior  of 
B  and  S  in  the  same  GB  (Tang,  Freeman  and  Olson  [6])  is  attached  as  Appendix  B. 

A  major  goal  of  our  research  is  the  exploitation  of  third-element  interactions  to  enhance 
GB  cohesion  and  minimize  embrittlement  susceptibility.  While  the  highly  precise  FLAPW  and 


DMol  calculations  require  use  of  CRAY  supercomputing  facilities  provided  by  NSF,  a  survey  of 
third-element  interaction  effects  is  being  conducted  employing  the  cluster  Discrete  Variational 
Method  (DVM)  which  can  be  run  on  the  HP730  Workstation  obtained  last  year  under  the  subject 
ONR  grant.  A  brief  description  of  its  application  to  P  embrittlement  is  included  in  a  recent 
review  paper  (Ellis  [7])  and  a  full  account  is  given  in  a  manuscript  (Sagert,  Ellis  and  Olson  [8]) 
provided  with  the  previous  report  which  explores  the  effect  of  Mo  on  the  P  embrittlement  of  the 
same  FeSSfllO]  (111)  GB. 

Our  recent  results  on  the  gettering  of  impurities  by  novel  compounds  in  rapidly  solidified 
steels  were  summarized  in  a  conference  abstract  paper  (Ghosh  and  Olson  [9])  also  provided  with 
the  previous  report. 

2.  EMBRITTLEMENT  ENERGETICS 

Table  1  of  Appendix  A  gives  a  comprehensive  summary  of  our  most  precise  FLAPW 
calculations  of  the  energy  contributions  to  the  effects  of  both  P  and  B  on  the  cohesion  of  the 
Fe  Z3  [iTO]  (111)  GB.  The  values  for  Fe-P  are  from  the  large  supercell  calculations  of  Ref.  [4] 
employing  the  GB  relaxations  of  Ref.  [3],  while  the  values  for  Fe-B  are  from  very  recently 
completed  calculations  employing  the  same  large  supercells  with  the  GB  relaxations  of  App.  B 
[6].  Values  are  listed  for  the  binding  energies  per  impurity  atom  (relative  to  an  isolated  impurity 
monolayer)  for  the  GB  (AEb  )  and  FS  (AEs)  for  (a)  unrelaxed  configurations  defining  the 
"chemical"  interaction,  (b)  the  energy  difference  associated  with  boundary  relaxation  defining 
the  "mechanical"  energy  penalty  of  boundary  reconstmction,  and  (c)  the  total  energy  change  for 
relaxed  configurations  taking  into  account  both  contributions.  Also  listed  are  the  values  for  non¬ 
spin-polarized  or  nonmagnetic  (NM)  calculations  to  compare  with  the  spin-polarized  or 
ferromagnetic  (FM)  calculations.  These  ultimately  define  the  magnetic  contribution  to  the 
AEb-AEs  energy  differences  governing  embrittlement  potency. 

The  Fe-P  results  for  the  unrelaxed  case  are  similar  to  the  earlier  small  cell  results  [1]  in 
that  the  pure  chemical  contribution  to  AEb-AEg  is  cohesion  enhancing  by  -0.75eV.  This  is  more 
than  offset,  however,  by  the  mechanical  reconstruction  energy  penalty  which  is  much  larger  in 
the  GB,  making  the  final  AEb-AEg  positive  and  therefore  embrittling  by  +1.60eV  in  the 
nonmagnetic  case  and  +0.79eV  in  the  ferromagnetic  case.  The  results  for  Fe-B  show  that 
compared  to  P  the  B  has  a  similar  binding  energy  in  the  GB  (with  or  without  relaxation)  but  a 
significantly  smaller  binding  energy  on  the  FS  consistent  with  the  earlier  calculations  of  Ref  [2]. 
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Experimental  estimates  of  AEb-AEs  are  +0.5eV  for  P  and  -0.8eV  for  B  (with  an 
uncertainty  of  -O.leV),  consistent  with  observed  behavior  of  P  as  an  embrittler  and  B  as  a 
cohesion  enhancer.  The  calculated  difference  between  P  and  B  is  reasonably  in  line  with 
experiment  but  the  individual  AEb-AEg  values  are  shifted  towards  embrittlement  by  0.3  to  0.8eV. 
This  may  be  due  to  a  systematic  error  in  the  calculated  relaxation  energies  which  are  the  most 
difficult  quantity  to  determine. 

Toward  improved  calculation  of  GB  relaxations  with  and  without  P  and  B,  revised  DMol 
cluster  calculations  have  now  been  initiated  with  modified  boundary  conditions.  The  cluster 
models  are  relaxed  at  constant  volume,  locking  cluster  surface  atoms  in  place  to  minimize  free 
surface  effects  arising  from  the  finite  cluster  size.  While  the  previous  calculations  employed  a 
volume  based  on  the  initial  CSL  GB  configuration,  the  new  calculations  employ  an  initial 
volume  taken  from  the  embedded-atom  method  (EAM)  calculations  of  Krasko  [10]  for  the  same 
boundaries.  This  allows  for  GB  free  volume  that  was  not  allowed  in  the  earlier  calculations,  and 
may  have  provided  an  effective  pressure  which  raised  the  relative  GB  energy  in  the  presence  of 
the  segregants.  The  new  calculations  actually  give  a  higher  cluster  energy  with  the  free  volume 
allowed,  but  we  believe  this  is  an  artifact  of  the  finite  cluster  size.  FLAPW  calculations  with  the 
new  relaxed  configurations  will  determine  if  the  GB  energy  is  now  lower  in  the  bulk 
environment.  We  are  hopeful  that  this  modification  will  resolve  the  apparent  systematic  error  in 
GB  relaxation  energies. 

While  calibration  of  our  calculations  against  experimental  energies  is  important  to  assess 
whether  our  methods  realistically  incorporate  the  essential  physics  of  the  problem,  most 
important  is  an  understanding  of  the  trends  between  species  from  which  insights  can  be  gleaned 
for  grain  boundary  design.  Fortunately,  the  chemical  bonding  trends  we  observe  are  much  less 
sensitive  to  cell  size  and  details  of  structural  relaxation. 

3.  CHEMICAL  BONDING 

To  assess  the  chemical  bonding  behaviors.  Figure  1  shows  computed  charge  transfer  plots 
comparing  P  and  B  in  the  GB  and  FS  environments.  These  are  obtained  by  subtracting  the 
superimposed  charge  densities  of  the  pure  components  from  the  computed  total  charge  density. 
Contours  start  from  1x10'^  e/a.u.^  and  increase  successively  by  a  factor  of  green,  yellow, 
red  and  pink  colors  denote  positive  differences,  and  light  and  dark  blue  denote  negative 
differences.  While  the  results  for  P  are  similar  to  those  obtained  in  our  earlier  small  cell 
calculations  [1]  the  plots  in  Figure  1  are  taken  from  the  large  supercell  calculations  of  Appendix 
A  [5]. 
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Figure  1 .  Calculated  valence  charge  density  difference  for  P  and  B  in  grain  bound^  (GB)  and 
free  surface  (FS)  environments.  Contours  start  from  1x10'^  e/a.u.^  and  increase 
successively  by  a  factor  of  2^^^ ;  green,  yellow,  red  and  pink  colors  denote  positive 
differences,  light  and  dark  blue  denote  negative  differences. 


A  deeper  interpretation  of  the  charge  transfer  plots  has  emerged  from  detailed  analysis  of 
specific  molecular  orbital  wave  functions  in  the  DMol  cluster  calculations  for  P  in  Ref.  [3]  and 
for  B  and  S  in  App.  B  [6].  As  illustrated  for  an  S  3p  orbital  (1^2)  in  Fig.  3a  of  App.  B  and  for  a 
B  2p  orbital  (Aii)  in  Fig.  6  of  App.  B,  we  find  that  the  3p  electrons  of  the  embrittlers  P  and  S 
actually  show  very  little  hybridization  with  Fe,  while  the  2p  electrons  of  B  do  hybridize  with  Fe 
to  give  significant  covalent  a  bonding  normal  to  the  GB  plane.  Although  the  local  charge 
transfer  maxima  for  P  in  Figure  1  resemble  covalent  bonds,  we  interpret  the  behavior  of  P  as  an 
"embedded"  atom  which  simply  interacts  electrostatically  with  the  potential  of  the  surrounding 
Fe  atoms,  the  local  charge  transfer  maxima  corresponding  predominantly  to  P3p  electrons  which 
are  not  significantly  hybridized  with  the  neighboring  Fe.  This  electrostatic  interaction  is 
relatively  unchanged  in  character  on  going  from  the  GB  to  the  FS  (although  the  P  3p  energy 
levels  move  in  response  to  the  changed  potential,  and  this  slightly  favors  interaction  in  the  GB 
environment).  In  contrast,  the  local  charge  transfer  maxima  for  B  do  reflect  directional  covalent 
cy  bonding,  an  important  consequence  of  which  is  the  "dangling  bond"  above  the  B  atom  on  the 
FS  as  revealed  in  Figure  1.  We  thus  associate  this  directional  bonding  character  of  B  with  the 
larger  chemical  contribution  to  AEb-AEs,  the  primary  contribution  to  its  differing  behavior 
relative  to  P. 

The  corresponding  spin  density  difference  plots  for  the  configurations  of  Figure  1  are 
shown  in  Figure  2,  employing  the  same  units  and  color  code.  Similar  to  the  earlier  small  cell 
calculations  of  Ref.  [1]  it  is  evident  that  the  impurities  reduce  the  magnetic  moment  (blue)  of 
near  neighbor  Fe  atoms  in  both  the  GB  and  FS  environments.  An  important  feature  revealed  by 
the  new  large  cell  calculations  is  that  Fe  atoms  further  from  the  impurities  have  enhanced 
moments  (red),  and  the  magnetic  effects  have  a  significantly  longer  range  than  the  net  charge 
transfer  effects.  The  change  of  magnetic  moments  of  the  atoms  of  figure  2  are  summarized  in 
Table  2  of  Appendix  A  [5].  Also  listed  are  the  moments  in  the  relaxed  clean  Fe  boundary  where 
the  CO  phase  transition  causes  an  antiferromagnetic  state  [3].  The  net  effect  of  the  longer  range 
magnetic  interactions  and  the  co  relaxation  of  the  clean  boundary  is  summarized  by  the  magnetic 
energy  increments  in  the  last  column  of  Table  1  of  Appendix  A  [5].  In  contrast  to  the  earlier 
small  cell  estimates  [1],  the  contributions  of  magnetism  actually  decrease  the  embrittlement 
potency  of  P  by  -O.SleV.  The  magnetic  effect  for  B  is  the  opposite,  decreasing  its  cohesion 
enhancing  effect  by  40.38eV. 
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Figure  2.  Calculated  spin  density  difference  for  P  and  B  in  grain  boundary  (GB)  and  free 
surface  (FS)  environments.  Units  and  color  code  as  in  Figure  1. 


Table  1 


Integrated  Charge  and  Energy  Densities  Between  Atoms 
(Values  taken  for  bond  radii  =^0.53  A, 

Excluded  length  =  0.26  A) 


Clean  Mo/Fe  Grain  Boundary  Clust^ 


Atoms) 


Vac(I)-Mo(ll) 

Vac(I)-Fe(lIl) 

Vac(I)-Fe(IV) 

Mo(II)-Mo(lI) 

Fe(Ill)-Fe(ni) 


Length 

1.81 

2.17 

1.82 

3.51 

2.17 


Total  Charj 
0.40 
0.89 
0.83 
1.00 
1.34 


P/Mn/Fe  Grain  Boundary  Cluster — (35 — Atoms)^ 


R  n  n  d 

0 

Bond  Length  (A) 

Total  Charge  (e) 

Total  Energy 

P(I)-Mo(Il) 
P(I)-Fe(III) 
P(I)-Fe(lV) 
Mo(l  I)-Mo(ll) 
Fe(llI)-Fe(lII) 

1.81 

2.17 

1.82 

3.51 

2.17 

0.80 

1.19 

1.09 

1  34 

1 .33 

-1.41 

-3.77 

-3.00 

-3.45 

-4.05 

Vac(I)-Mo(lI) 
Vac(I)-Fe(lII) 
V  ac(I)-Fe(I  V) 
Mo(II)-Mo(II) 


P/Mo/Fe 


Free  Surface  Cluster  (36  Atoms) 


Bond 

Bond  Length  (A) 

Total  Charge  (e) 

Total  Energ 

P(I)-Mo(lI) 

P(I)-Fe(lIl) 

P(I)-Fe(lV) 

Mo(II)-Mo(ll) 

1.86 

2.15 

1.61 

3.52 

0.98 

1.20 

1.10 

1.50 

-0.36 

-2.16 

-0.94 

-1.59 

4.  ALLOYING  EFFECTS 


As  described  in  the  previous  progress  report,  the  DVM  cluster  method  is  being  applied  to 
the  third-element  interaction  with  Mo  employing  the  same  atomic  configurations  as  in  Fig.  1,  but 
with  a  substitutional  monolayer  of  Mo  at  the  GB  core.  The  Mo-P  interaction  in  the  intergranular 
embrittlement  of  Fe  is  fairly  well  studied  experimentally,  the  evidence  indicating  that  grain 
boundary  cosegregation  of  Mo  and  P  results  in  a  reduction  of  embrittlement.  Although  the 
charge  transfer  behavior  looks  quite  different  from  the  binary  Fe-P  case  of  Figure  1 ,  the  P  is 
found  to  still  show  "embedded"  behavior  with  no  significant  hybridization  with  either  Fe  or  Mo. 
There  is  however  enhanced  charge  transfer  from  Mo  to  P  in  association  with  its  greater 
electronegativity  difference,  resulting  in  a  greater  ionic  bonding  component  and  an  electrostatic 
attraction  of  the  P  3p  electrons  to  the  region  adjacent  to  the  Mo.  The  behavior  is  similar  in  the 
GB  and  FS  environments,  but  with  slightly  greater  charge  transfer  (and  ionic  interaction)  on  the 
FS.  The  energetic  result  is  that  Mo  attracts  P  in  both  environments,  but  slightly  more  on  the  FS 
giving  an  increased  chemical  contribution  to  AEb-AEs  for  P  of  +0.8eV.  However,  this  apparent 
increased  embrittlement  potency  of  P  in  the  presence  of  Mo  is  found  to  be  more  than 
compensated  by  the  direct  cohesion  enhancing  effect  of  Mo.  Hence  the  attraction  of  Mo  to  the 
GB  by  segregated  P  (with  a  Mo-P  interaction  energy  of  ~-leV  per  atom)  will  overall  enhance 
boundary  cohesion. 

The  origin  of  the  Mo  cohesion  enhancement  is  being  investigated  in  terms  of  detailed 
local  bonding.  Labelling  atomic  positions  according  to  the  scheme  of  Figure  3  with  Mo 
occupying  site  n.  Table  1  summarizes  integrated  bond  charges  and  energies.  The  bond  volume 
is  defined  by  a  cylinder  of  0.53A  radius  with  a  length  determined  by  excluding  an  0.26A  distance 
from  the  atomic  cores.  Bond  energy  is  seen  to  correlate  well  with  bond  charge  with  the  notable 
exception  of  the  Mo  (n)-Mo(II)  bonds  which  have  an  anomolously  small  bond  energy.  This  is 
particularly  noticeable  for  the  clean  Mo/Fe  free  surface  cluster  where  the  Mo(n)-Mo(II)  bond  has 
the  highest  bond  charge  (1.29e)  but  second  smallest  bond  energy  of  -0.66eV.  Much  of  the 
charge  associated  with  these  bonds  apparently  involves  unfavorable  energetic  states  relative  to 
the  Mo  ground  state,  giving  in-plane  metal-metal  bonds  which  are  less  energetically  stable  than 
the  original  Fe(n)-Fe(II)  bonds.  That  this  behavior  is  more  pronounced  on  the  FS  than  the  GB 
may  be  very  significant  to  GB  cohesion.  Unusual  behavior  of  Mo  in  the  GB  environment  was 
also  reported  by  Dr.  G.  Krasko  of  ARL  at  our  SRG  subgroup  meeting  at  Harvard  in  December, 
based  on  LMTO-ASA  calculations.  This  behavior  will  be  investigated  further  and  compared 
with  that  of  embrittling  Mn  to  more  clearly  establish  the  electronic  features  of  cohesion 
enhancing  alloying  elements  for  alloy  design. 
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5.  EXPERIMENTS 


Results  described  in  our  previous  progress  report  demonstrated  that  with  careful 
avoidance  of  Si  contamination,  the  rapid  solidification  of  a  model  0.4C  Ni-Mo  steel  after  late 
addition  of  La  to  a  melt  deoxidized  with  Zr+Ti  successfully  reproduced  our  earlier  finding  of  La- 
P-O  compounds  which  efficiently  getter  P.  Rather  than  the  highly  stable  La202S  found  earlier, 
however,  the  S  was  found  to  be  gettered  in  the  form  of  TiS.  Analytical  microscopy  identified  the 
P  gettering  phases  as  predominantly  LaP04  and  LavPsOig-  Further  study  of  material  deoxidized 
with  Zr+Al  revealed  La202S  as  the  dominant  S  gettering  phase,  and  led  to  the  discovery  of  ZrP 
as  a  P  gettering  phase  as  reported  in  Ref.  [9].  Further  correlation  of  gettering  phases  with  melt 
deoxidation  practice,  aided  by  thermodynamic  modeling,  will  address  predictive  control  of 
impurity  gettering. 

As  part  of  his  surface  analytical  study  of  embrittling  impurities  under  AASERT  support. 
Graduate  Student  David  Spaulding  has  measured  the  grain  boundary  fracture  strength  in 
polycrystalline  tensile  specimens  of  model  UHS  steels  in  which  he  previously  measured  the  grain 
boundary  segregation  of  P  and  S  as  a  function  of  tempering  time  by  Scanning  Auger 
Microanalysis.  The  high  yield  strength  of  these  alloys  allowed  fracture  strength  to  be  measured 
under  well  defined  homogeneous  elastic  loading  conditions.  Correlation  of  boundary  strength 
with  the  predicted  reduction  in  ideal  thermodynamic  work  of  fracture  will  allow  a  more  precise 
test  of  the  Rice-Wang  model  than  previous  correlations  with  shifts  in  ductile-brittle  transition 
temperature.  Further  experiments  will  quantify  the  effect  of  hydrogen  charging. 
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The  effects  of  P  and  B  impurities  on  the  thermodynamic,  electronic  and  mag¬ 
netic  properties  of  both  the  free  Fe(lll)  surface  and  FeE3[ll0]  (111)  gram  boundary 
(GB)  are  investigated  using  the  local  density  full  potential  linearized  augmented  plane 
wave  method.  We  find  that  the  spatial  anisotropy  of  the  impunty-Fe  bonding  plays  a 
key  role  in  the  relative  cohesion  of  the  FeE3[ll0]  (111)  GB.  The  segregation  energies 
of  P  and  B  (relative  to  an  isolated  monolayer)  are  almost  equal  in  the  GB  environ¬ 
ment  even  when  the  impurity-induced  structural  relaxation  in  the  FeE3[nO]  (111) 
GB  is  considered.  For  P,  the  segregation  energy  at  the  free  surface  is  lower  than  at 
the  GB  by  -0.81  eV  mainly  due  to  the  large  relaxation  energy  in  Fe  GB.  This  energy 
difference  is  consistent  with  the  observed  embrittlement  potency  of  P  as  described 
by  the  Rice- Wang  thermodynamic  model  of  Gb  embrittlement.  The  bonding  of  P  in 
both  environments  is  electrostatic  in  nature,  shows  no  significant  hybridization  with 
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Fe.  In  contrast,  B  shows  strong  hybridization  with  Fe  giving  directional  bonding  nor¬ 
mal  to  the  GB  plane  which  increases  the  relative  energy  on  the  Fe(lll)  surface,  m 
line  with  the  relative  cohesion  enhancing  effect  of  B.  Magnetic  contributions  reduce 
both  the  embrittling  effects  of  P  and  the  cohesion  enhancing  effect  of  B. 
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I.  INTRODUCTION 


A  thermodynamic  theory  developed  by  Rice  and  Wang  [1,2]  describes  the  mech¬ 
anism  of  metalloid-induced  intergranular  embrittlement  through  the  competition  be¬ 
tween  dislocation  crack  blunting  and  brittle  boundary  separation.  While  crystal 
plasticity  considerations  show  interesting  directional  effects  on  the  relative  ease  of 
dislocation  emission  from  the  crack  tip,  tested  in  critical  bicrystal  experiments,  the 
most  striking  result  of  the  analysis  is  the  prediction  that  the  potency  of  a  segregat¬ 
ing  solute  in  reducing  the  “Griffith  work”  of  brittle  boundary  separation  is  a  linear 
function  of  the  difference  in  segregation  energy  for  that  solute  at  a  grain  boundary 
(GB)  AEbi  and  a  free  surface  (FS),  AEs-  Simply  put,  a  solute  with  a  positive  energy 
difference  AEt,  -  AEs  (i.e.,  AEb  is  less  negative)  will  be  a  more  potent  embrittler, 
or  vice  versa.  The  validity  of  this  hypothesis  has  been  supported  experimentally  for 
Fe-base  materials  for  which  significant  of  data  exist  for  both  the  relative  embrittling 
potencies  of  solutes  and  the  relevant  surface  thermodynamic  quantities  [3]. 

A  predictive  theory  of  intergranular  cohesion  can  be  developed  by  directly  deter¬ 
mining  both  AEb  and  AE^  on  an  equal  footing  with  present  state-of-the-art  energy 
band  approaches;  this  offers  the  promise  of  major  improvements  in  boundary  cohe¬ 
sion  sensitive  properties  such  as  the  stress  corrosion  resistance  of  ultra-high  strength 
steels.  To  this  end,  we  have  carried  out  full-potential  linearized  augmented  plane  wave 
(FLAPW)  total  energy  calculations  to  investigate  the  effects  of  a  phosphorus  impu¬ 
rity  in  the  FeS3[ll0]  (111)  GB  and  on  the  corresponding  Fe(lll)  FS  [4].  In  earlier 
calculations  a  cell  that  contains  only  the  first  neighbor  Fe  atoms  to  the  P  impurity, 
we  found  that  the  Fe-P  chemical  interaction  is  slightly  stronger  at  the  Fe(lll)  FS  due 
to  a  shorter  interatomic  distance.  However,  as  a  result  of  the  small  size  of  the  unit 
cell  and  the  exclusion  of  impurity-induced  structural  relaxation,  the  calculated  AEb  - 
AEs  was  found  to  be  negative  (-1.1  eV/cell),  in  contrast  with  the  known  embrittling 
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potency  of  P  in  Fe. 

In  a  subsequent  calculation,  we  expanded  the  unit  cell  size  for  more  realistic 
simulations  (24  atoms/cell  for  the  GB  and  15  atoms/cell  for  the  FS)  [14].  By  using 
the  equilibrium  geometry  determined  from  DMol  cluster  calculations  [8],  we  found 
that  structural  relaxation  may  reduce  the  total  energy  of  the  clean  FeS3[ll0]  (111) 
(ie.,  the  reference  system  for  AE,)  by  about  1.61  eV/cell.  This,  in  turn,  results  in  a 
large  correction  to  AEb-  By  contrast,  the  impurity  atoms  do  not  induce  significant 
reconstruction  for  the  less  constrained  environment  of  the  Fe(lll)  surface  [5]  and  thus 
AEs  is  only  slightly  affected.  As  a  result,  AEb  -  AEs  becomes  positive,  0.79  eV  per 
P  atom  -  in  good  agreement  with  experimental  data. 

In  the  present  paper,  we  explore  the  mechanism  of  impurity-induced  embrittle¬ 
ment  more  deeply  by  comparing  the  behaviors  of  P  and  B  in  the  same  FeE3[ll0]  (111) 
GB  and  Fe(lll)  FS  environments.  As  B  is  known  to  be  a  GB  cohesion  enhancer  m 
Fe  [1-3],  the  comparison  offers  insights  into  the  key  electronic  features  underlying  the 
influence  of  segregants  on  relative  boundary  cohesion.  Following  presentation  of  the 
model  and  computational  details  in  Sec.  II,  we  discuss  the  difference  between  B  and 
P  on  the  Fe  boundary  cohesion  from  the  point  of  view  of  the  calculated  electronic 
and  magnetic  properties  in  Sec.  III. 

II.  MODEL  AND  COMPUTATION 

As  sketched  in  Fig.  1,  we  used  a  slab  model  to  simulate  both  the  P/Fe(lll) 
free  surface  and  P/FeE3[ll0]  (111)  GB,  which  minimizes  GB  interactions  inherent  in 
superlattice  cells  [4,6].  For  the  B/Fe(lll)  and  P/Fe(lll)  FS  systems,  the  Fe(lll) 
substrate  is  simulated  by  a  13  layer  slab  and  the  B  and  P  adsorbates  are  placed 
pseudomorphically  on  the  three-fold  hollow  sites  on  both  sides  of  the  slab.  For  the 
GB  systems,  a  23  layer  slab  is  adopted  to  simulate  the  clean  FeE3[lT0]  (111).  With  12 
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FIG.  1.  Model  and  notation  for  the  structure  of  Fe  and  impnrities  (I)  at  (a)  Fe 

i;3[ll0]  (111)  “incoherent”  twin  boundary,  and  (b)  Fe(lll)  free  surface. 


Fe  layers  in  between,  the  interaction  between  the  free  surface  (introduced  artificially 
in  the  slab  model)  and  the  GB  is  expected  to  be  sufficiently  reduced.  The  2D  lattice 
constant  and  the  unrelaxed  Fe-Fe  interatomic  distance  are  chosen  from  experimental 
values  for  bulk  bcc  Fe,  i.e.,  dfe-Fe=4.69  a.u.  The  slab  models  are  sufficiently  thick 
to  allow  multilayer  atomic  relaxation  for  Fe  in  both  FS  and  GB  environments.  For 
the  clean  and  impurity  adsorbed  Fe(lll)  surface,  the  equilibrium  atomic  geometry 
is  determined  through  FLAPW  total  energy  minimization  [5,7].  For  gram  boundary 
systems,  geometries  obtained  from  DMol  cluster  force  calculations  [8]  are  adopted. 

In  the  FLAPW  method  [9],  no  shape  approximations  are  made  to  the  charge 
densities,  potentials  and  matrix  elements.  The  core  states  are  treated  fully  relativis- 
tically  and  the  valence  states  are  treated  semi-relativistically  (i.e.,  without  spin-orbit 
coupling)  [10].  We  employ  the  Hedin  and  Lundqvist  and  the  von  Barth-Hedin  formu¬ 
las  for  the  exchange-correlation  potentials  for  the  nonmagnetic  and  the  spin-polarized 
calculations,  respectively  [11].  This  approach  has  been  applied  very  successfully  in  the 
last  decade  to  determine  the  electronic  and  magnetic  properties  of  many  transition 
metal  systems.  [12] 

Energy  cutoffs  of  11  Ry  and  70  Ry  are  employed  for  plane  wave  bases  and  star 
functions  to  describe  the  wave  functions  and  the  charge  density  and  potential  in  the 
interstitial  region,  respectively.  Within  the  muffin-tin  (MT)  spheres  {rMT,Fe  =  2.0 
a.u.,  T'MT p  ~  I'S  a.u.,  tmt,b  ~  a.u.),  lattice  harmonics  with  angular-momentum 

/  up  to  8  are  adopted.  Convergence  is  assumed  when  the  average  root-mean-square 
difference  between  the  input  and  output  charge  (spin)  densities  is  less  than  5  x  10  ^ 
e/(a.u.)^  (1  X  10“^  e/(a.u.)^).  The  step-forward  fixed-basis  approach  [13]  is  used  to 
speed  up  the  calculations. 


III.  CALCULATED  RESULTS 
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A.  Structural  relaxation 


In  previous  calculations  [14,8],  we  found  that  the  clean  FeE3[ll0]  (111)  GB 
undergoes  an  u)  phase  transition  which,  in  turn,  results  in  an  antiferromagnetic  cou¬ 
pling  in  the  core  region  of  the  GB.  In  addition,  the  P  impurities  strongly  push  the 
Fe  atoms  apart  across  the  GB  to  form  an  “anti-u;”  structure,  corresponding  to  a 
metastable  state  of  the  clean  FeS3[ll0]  (HI)  GB.  As  listed  in  Table  1,  such  a  large 
reconstruction  gives  a  large  “mechanical  energy”  contribution  to  the  final  value  of 
by  1.61  eV/cell.  The  mechanical  energy  contribution  corresponds  to  the  energy 
of  the  unrelaxed  Fe  boundary  after  impurity  removed  ,  relative  to  the  relaxed  clean 
Fe  boundary. 

Although  the  B  impurity  is  significantly  smaller  than  the  P  atom,  it  still  cannot 
be  well  accommondated  in  the  a;-Fe  GB  and  some  B-induced  Fe  reconstruction  should 
also  be  expected.  Indeed,  the  calculated  B-Fe(3)  bond  length  is  4.06  a.u.,  only  9% 
shorter  than  that  between  P-Fe(3),  4.42  a.u.  [8].  Surprisingly,  this  smaller  relaxation 
gives  rise  to  an  even  larger  mechanical  energy  in  Table  1,  1.69  eV,  compared  to  that 
in  the  P /Fe  case  for  which  the  Fe  structure  becomes  closer  to  the  metastable  anti-w 
structure. 

As  found  previously  [4],  the  bond  length  is  shorter  in  the  FS  environment  than 
that  in  the  GB;  the  calculated  dp-Fe(z)  iii  P/Fe(lll)  is  4.5  a.u.  and  dp-Feiz)  in 
B/Fe(lll)  is  3.44  a.u.  However,  since  the  impurity  atoiiis  have  a  large  space  for 
adjustment  in  the  vacuum  region,  the  mechanical  effects-are  expected  to  be  less  im¬ 
portant  in  the  FS  environment.  Indeed,  as  also  listed  in  Table  1,  the  mechanical 
energy  for  the  surface,  is  usually  less  than  0.2  eV/cell,  even  when  the  contribu¬ 
tions  from  both  the  impurity  relaxation  and  the  impurity-induced  Fe  relaxation  are 
considered  [14]. 
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TABLE  I.  Calculated  binding  energies  (in  eV)  of  P  and  B  in  the  Fe  GB  and  FS  envi 
ronment  and  the  decompositions  into  chemical,  mechanical  and  magnetic  contributions 


AEb 

AEs 

AEf,  -  AEs 

Magnetic 

P-Fe 

Unrelaxed  (Chemical) 

NM 

-8.16 

—  7.58 

-0.58 

-0.17 

FM 

-8.19 

-7.44 

-0.75 

Relaxation  (Mechnical) 

NM 

+2.19 

+0.01 

+2.18 

-0.64 

FM 

+  1.61 

+0.07 

+1.54 

Relaxed  (Total) 

NM 

-5.97 

-7.57 

+  1.60 

-0.81 

FM 

-6.58 

-7.37 

+0.79 

B-Fe 

Unrelaxed  (Chemical) 

NM 

-8.27 

-6.46 

-1.81 

+0.30 

FM 

-8.18 

-6.67 

-1.51 

Relaxation  (Mechnical) 

NM 

+  1.53 

+0.10 

+1.43 

+0.08 

FM 

+  1.69 

+0.18 

+1.51 

Relaxed  (Total) 

NM 

-6.74 

-6.36 

-0.38 

+0.38 

FM 


-6.49 


-6.49 


+0.00 


B.  Chemical  interaction 


The  impurity  induced  charge  redistributions  for  the  P/Fe(lll)  and  B/Fe(lll) 
FS  systems,  obtained  by  subtracting  the  superimposed  charge  density  from  a  free  P 
or  B  monolayer  and  a  clean  Fe(lll)  surface  from  the  self-consistent  charge  density  for 
the  corresponding  adsorption  system,  are  presented  in  Figs.  2(a)  and  2(b),  respec¬ 
tively.  Contours  with  green,  yellow,  red  and  pink  colors  denote  the  positive  charge 
density  difference  and  those  with  light  and  dark  blue  colors  represent  the  negative 
charge  density  difference.  For  P/Fe(lll),  strong  charge  accumulations  are  shown  in 
the  region  between  both  P-Fe(3)  and  P-Fe(l)  -  indicating  some  chemical  interaction 
between  these  atoms.  The  charge  density  at  the  inner  region  around  the  P  atom  is 
significantly  decreased  due  to  the  effects  of  the  surrounding  Fe  atoms.  This  appar¬ 
ent  reverse-charge  transfer  contradicts  a  simple  estimate  from  the  electronegativity 
(2.19  for  P  and  1.83  for  Fe).  As  demonstrated  previously,  these  behaviors  can  be 
understood  from  the  large  spatial  extension  of  the  P-3p  wave  function,  and  thus  the 

“embedding”  character  of  the  P-Fe  bonding  [14]. 

By  contrast,  since  the  B-2p  wave  function  is  more  localized  than  the  P-3p  wave 
function,  the  B-Fe(3)  bonding  is  more  covalent-like.  In  addition,  B  has  only  one  2p 
electron  while  P  possesses  three  3p  electrons.  In  the  isolated  monolayer  geometry, 
the  B-2p  electron  occupies  the  bonding  p,;,y  states  which  leaves  the  non-bonding 
state  empty.  However,  the  pz  state  becomes  lower  in  energy  by  hybridizing  with  the 
Fe(3)-d^2  state  and  thus,  as  shown  in  Fig.  2(b)  there  is  charge  transfer  from  the  in¬ 
plane  px,y  states  to  the  vertical  pz  state.  Therefore,  B-Fe  bonding  in  Fig.  2(b)  shows 
much  stronger  spatial  anisotropy  compared  to  that  between  P-Fe  in  Fig.  2(a),  i.e., 
stronger  vertical  B-Fe(3)  bonding  and  weaker  lateral  B-Fe(l)  bonding.  As  a  result, 
the  chemical  part  of  for  B  in  the  FS  environment,  6.67  eV/adatom,  is  0.77  eV 
smaller  than  that  for  P. 
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(a)  (b) 


FIG.  2.  The  calctilated  valence  charge  density  difference  for  the  (a)  P/Fe(lll)  and 
(b)  B/Fe(lll)  FS  systems.  Contours  start  from  1  x  10"^  e/a.u.^  and  increase  successively 
by  a  factor  of  2« ;  green,  yellow,  red  and  pink  colors  denote  positive  differences  and  light 
and  dark  blue  colors  denote  negative  differences. 


The  correspopding  charge  density  differences  for  P/Fe  and  B/Fe  GB  systems 
are  presented  in  Figs.  3(a)  and  3(b),  respectively.  Due  to  the  longer  bond  length, 
the  charge  accumulation  and  hence  the  strength  of  Fe-impurity  bonding  are  weaker 
compared  to  the  case  in  the  FS  environment.  However,  there  are  two  P(B)-Fe(3) 
bonds  across  the  GB,  which  results  in  a  larger  chemical  energy  in  the  GB  systems 
by  0.75  eV  for  P  and  1.51  eV  for  B.  Dramatically,  although  the  B-Fe(l)  interaction 
is  much  weaker  compared  to  P-Fe(l)  for  the  reason  discussed  above,  the  chemical 
energies  for  B  and  P  in  the  GB  environment  become  almost  ec[ual,  8.18  eV. 

As  found  previously  [4, 15],  the  presence  of  the  impurity  reduces  the  Fe(2)-Fe(2) 
interaction  across  the  GB.  In  Fig.  2,  the  negative  charge  density  difference  can  also 
be  found  above  the  Fe  atoms  in  the  FS  environment.  This  charge  removal  was  once 
proposed  to  explain  the  embrittlement  effect  of  impurities  on  GB  cohesion  [15].  How¬ 
ever,  since  this  effect  is  even  stronger  for  B  which  is  a  known  GB  cohesion  enhancer 
in  Fe,  this  does  not  appear  to  be  an  essential  feature  of  embittlement  behavior. 

C.  Magnetic  interaction 

In  the  previous  small  cell  calculations  [4],  we  found  that  P  impurities  strongly 
diminish  the  magnetization  of  their  surrounding  Fe  atoms.  This  is  still  the  case 
in  the  present  large  cell  calculation.  As  shown  in  Fig.  4(a)  by  the  spin  density 
difference  contours  for  the  P/Fe  GB,  the  detrimental  effects  of  the  P  impurities  on 
the  magnetization  around  the  Fe(l),  Fe(2)  and  Fe(3)  atoms  are  very  clear.  However, 
P  induces  an  enhancement  of  magnetization  around  the  second-rank  Fe  atoms,  that 
was  not  obtained  in  the  previous  small  cell  calculations  [4].  Note  that  this  behavior 
is  not  seen  in  Fig.  3(a)  for  the  charge  density  difference  —  indicating  that  the  P 
induces  a  short-ranged  perturbation  (screening)  for  the  charge  distribution  but  shows 
a  long-ranged  oscillatory  behavior  for  the  magnetic  disturbance.  The  spin  density  at 
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FIG.  3.  The  calculated  valence  chargi 
GB.  Contours  start  from  1  X  10"^  e/a.u.^  am 
yellow,  red  and  pink  colors  denote  positive  di 


negative  differences. 


(a)  (b) 


FIG.  4.  The  calculated  spin  density  difference  for  (a)  P/Fe  and  (b)  B/Fe  GB.  Con¬ 
tours  start  from  1  x  10“^  e/a.u.^  and  increase  successively  by  a  factor  of  2« ;  green,  yellow, 
red  and  pink  colors  denote  positive  differences  and  light  and  dark  blue  colors  denote  negative 


differences. 


the  P  site  is  very  small,  so  P  does  not  carry  a  large  magnetic  moment  (ss-0.01  fia) 
although  its  3p  band  exhibits  some  exchange  splitting  [14].  By  comparison,  as  shown 
in  Fig.  3(b),  B  effects  on  the  magnetization  of  the  surrounding  Fe  atoms  are  weaker, 
especially  for  the  Fe(2)  atom.  In  addition,  the  magnitude  of  the  negative  spin  density 
around  the  B  site  is  not  small.  B  carries  a  magnetic  moment  of  -0.05  ps  m  its 
muffin-tin  sphere.  In  Table  1,  we  can  see  from  comparison  of  ferromagnetic  (FM)  and 
nonmagnetic  (NM)  calculations  (i.e.,  with  and  without  spin  polarization)  that  the 
magnetic  contributions  to  the  segregation  energies  (magnetic  energy)  are  very  small, 
typically  around  0.1  eV.  However,  this  can  be  important  in  some  case  such  as  for  the 
mechanical  energy  in  the  P  case. 

The  spin  density  differences  plotted  in  Figs.  5(a)  and  5(b)  for  the  P/Fe(lll)  and 
B/Fe(lll)  FS  system  show  the  same  trend  as  that  in  Fig.  4,  i.e.,  P  and  B  adsorbates 
diminish  (enhance)  the  magnetization  around  their  first  (second)  rank  Fe  neighbors. 
However,  the  difference  in  magnitude  is  also  obvious  for  both  P  and  B.  For  P,  the 
spin  density  around  Fe(2)  and  Fe(l)  is  less  affected  in  the  FS  environment,  although 
the  interatomic  distances  are  smaller  compared  to  those  in  the  GB.  By  contrast,  B 
affects  Fe(2)  and  Fe(3)  more  strongly  in  the  FS.  As  reflected  in  Table  1,  the  magnetic 
energy  increment  (with  structural  relaxation)  for  P/Fe(lll)  is  0.20  eV  vs.  -0.13  eV 
for  B/Fe(lll),  while  it  is  -0.61  eV  for  the  P/Fe  GB  vs  0.25  eV  for  the  B/Fe  GB. 

Quantitatively,  the  calculated  magnetic  moment  in  each  Fe  muffin-tm  sphere 
is  presented  in  Table  2.  The  effects  of  P  are  obvious  since  it  reduces  the  magnetic 
moment  of  Fe(l)  by  1.0  ps  in  both  the  GB  and  FS  environments.  In  P/Fe(lll),  the 
P-induced  decrease  of  the  Fe(3)  magnetic  moment  is  also  as  large  as  0.63  pB,  which 
results  in  the  positive  magnetic  contribution  to  AE^.  By  comparison,  the  reduction  of 
the  Fe(3)  magnetic  moment  is  even  a  little  larger  in  B/Fe(lll),  0.65  hb-  By  contrast, 
B  effects  for  the  Fe(l)  magnetic  moment  are  reduced  in  the  FS,  0.38  pa  vs  0.45  pa 
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TABLE  II.  Calculated  magnetic  moments  (in  ms)  in  Fe.  GBO  represents  the  relaxe 
clean  Fe  GB,  while  GBl  and  GB2  denote  the  unrelaxed  GB  for  P  and  B,  respectively.  _ 


atom 

Fe  GBO 

Fe  GBl 

P/Fe  GB 

Fe  GB2 

Fe(l) 

-1.28 

2.40 

1.40 

2.50 

Fe(2) 

1.56 

2.11 

1.81 

2.18 

Fe(3) 

1.43 

1.81 

1.49 

2.07 

Fe(4) 

-1.75 

1.81 

2.01 

1.86 

Fe(5) 

0.80 

2.10 

2.05 

2.14 

B/FeGB  Fe(lll)  P/Fe(lll)  B/Fe(lll) 


2.05 

2.71 

1.74 

2.33 

2.17 

2.29 

2.22 

2.28 

1.90 

2.36 

1.73 

1.71 

2.15 

2.20 

2.30 

2.31 

2.14 

2.20 

2.17 

2.13 

Fe(6) 


1.42 


1.85 


1.98 


1.85 


2.02 


2.11 


2.27 


2.29 


(a)  (b) 


FIG.  5.  The  calculated  spin  density  difference  for  (a)  P/Fe(lll)  and  (b)  B/Fe(lll) 
FS.  Contours  start  from  1  x.  10“®  e/a.u.^  and  increase  successively  by  a  factor  of  2« ;  green, 
yellow,  red  and  pink  colors  denote  positive  differences  and  light  and  dark  blue  colors  denote 
negative  differences. 


in  the  GB;  this  contributes  a  negative  magnetic  energy  increment  for  of  B  since 
there  are  three  Fe(l)  atoms  but  only  one  Fe(3)  atom  at  the  surface.  However,  the 
magnetic  moment  enhances  in  both  B/Fe  GB  and  P /Fe  GB  if  compared  to  the  the 
strongly  depressed  value  for  the  relaxed  clean  Fe  GB  due  to  the  shorter  interatomic 
distance  and  the  antiferromagnetic  coupling  [14]. 

IV.  SEGREGATION  ENERGY  AND  DISCUSSION 

Finally,  the  calculated  results  for  the  AEh  -  AEs  difference  controlling  embrit¬ 
tlement  are  listed  in  Table  1.  Note  that,  by  including  the  impurity-induced  structural 
relaxation,  the  segregation  energy  for  P  in  the  grain  boundary  environment  is  smaller 
(less  negative)  than  that  on  the  free  surface,  i.e.,  AEi,  -  AEs  =  0.79  eV  >  0,  m 
good  agreement  with  experiment  (0.4±0.2  eV/cell).  According  to  the  Rice- Wang 
thermodynamic  theory  [1],  P  is  thus  an  embrittler  for  the  Fe£3[ll0]  (111)  GB.  By 
contrast,  AEb  -  AEs  is  almost  zero  (negative  for  the  nonmagnetic  case)  for  B,  and 
thus  B  is  weakly  cohesion  enhancing.  To  our  knowledge,  this  is  the  first  quantitative 
theoretical  explanation  for  the  different  behavior  of  B  and  P  metalloid  impurities  on 
GB  cohesion  in  Fe. 

From  the  discussions  above,  many  factors  affect  the  final  results  of  AEb  —  AE^- 
The  mechanical  energy  is  especially  important,  as  the  segregation  energy  differences 
are  negative  for  both  P  and  B  when  relaxation  for  the  clean  Fe  GB  is  neglected.  From 
Table  1,  we  can  see  that  the  magnetic  energies  are  not  large  for  either  AEb  or  AEs 
separately,  but  are  indeed  important  for  their  difference.  Therefore,  a  highly-precise 
spin  polarized  method  and,  importantly,  equal  treatment  for  the  FS  and  GB  systems, 
are  essential  to  obtain  correct  predictive  results.  The  net  magnetic  contributions  to 
the  AEb  —  AEs  differences  are  summarized  in  the  last  column  of  Table  1.  While  mag¬ 
netism  reduces  the  embrittlement  potency  of  P,  it  also  reduces  the  cohesion  enhancing 
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effect  of  B. 

Impurity-induced  reduction  of  Fe-Fe  bonding  is  here  eliminated  as  the  mech¬ 
anism  for  GB  embrittlement.  Since  the  embrittlement  behavior  of  an  impurity  is 
determined  by  AEb  -  it  is  essential  to  compare  the  difference  of  effects  in  the 

FS  and  GB  systems,  respectively.  In  both  the  FeS3[ll0]  (111)  GB  and  the  Fe(lll) 
FS,  there  are  three  Fe(l)  atoms.  The  strength  of  impurity-Fe(l)  bonding,  seen  from 
Figs.  2  and  3,  does  not  change  greatly  from  the  GB  to  FS  environments  and  thus 
is  not  significant  for  AEg  -  AEs-  By  contrast,  impurity-Fe(3)  bonding  is  expected 
to  play  the  important  role  since  one  out  of  two  of  the  vertical  bonds  is  broken  from 
the  GB  to  FS.  Therefore,  the  spatial  anisotropy  of  the  bonding  interaction  between 
the  impurity  states  with  the  surrounding  Fe  atom  can  be  crucial  to  embrittlement 
behavior.  Simply,  impurities  with  stronger  vertical  and  weaker  lateral  bonding  are 
favorable  to  enhance  the  GB  cohesion. 

This  conclusion  is  supported  by  the  comparison  of  P  and  B.  In  Fig.  3(b),  the 
B-Fe(3)  bonding  is  much  stronger  than  the  B-Fe(l)  bonding.  There  is  even  a  dangling 
bond  above  of  the  B  atom  for  B/Fe(lll)  in  Fig.  2(b).  Therefore,  B  prefers  the  GB 
environment  to  saturate  the  B-Fe(3)  bonding.  By  contrast,  P-Fe  bonding  is  more 
embedded-like  (electrostatic)  and  thus  there  is  no  such  dangling  bond  above  P  for 
P/Fe(lll)  in  Fig.  2(a).  In  addition,  the  strength  of  the  P-Fe(l)  bonding  is  almost  as 
strong  as  that  of  P-Fe(3)  bonding.  Only  one  out  of  five  P-Fe  bonds  is  broken  from 
the  GB  to  FS  for  P  (one  out  of  two  for  B)  and  thus  the  reduction  of  the  chemical 
energy  for  P  is  smaller  than  for  B.  This  is  the  main  reason  why  the  chemical  part 
of  AEb  —  AEs  for  B  is  two  times  larger  than  that  for  P.  Calculations  for  C  and  S 
impurities  are  in  progress  to  verify  the  generality  of  these  bonding  trends. 
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Abstract 


The  structure  and  electronic  properties  of  B  and  S  impurities  in  an  Fe  E3[110](lll)  grain  boundary 
are  studied  using  the  DMol  molecular  cluster  method.  A  large  cluster  containing  91  atoms  is  used 
to  simulate  the  local  environment  of  the  boundary.  Optimized  atomic  geometries  are  obtained  by 
atomic  force  calculations.  The  results  show  that  impurities  induce  large  relaxations  in  the  grain 
boundary,  and  that  the  electronic  structure  of  the  hybridized  B/Fe  system  is  very  different  from 
that  of  the  non-hybridized  S/Fe  system. 

PACS  Numbers:  61.70.-r,  61.70.Ng,  71.55.-i 
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I.  Introduction 


It  is  weU  known  that  impurities  such  as  P,  S,  Sn  and  Sb  cause  intergranular  embrittlement  of 
Fe,  while  C  and  B  enhance  the  intergranular  cohesion  in  iron  and  alloy  steels.'  Because  of  its  im¬ 
portance  to  grain-boundary  sensitive  properties  such  as  the  stress- corrosion  resistance  of  ultrahigh- 
strength  steels,  the  mechanism  underlying  these  phenomena  has  been  the  subject  of  considerable 
research  effort. It  has  become  evident  from  these  studies  that  the  desired  understanding  will 
require  detailed  consideration  of  precise  grain  boundary  atomic  structure  and  energetics  coupled 
with  underlying  electronic  level  structure  and  properties. 

Recent  theoretical  efforts  have  applied  first  principles  electronic  structure  calculations  to  study 
the  role  of  impurities  in  the  embrittlement  of  Fe  and  proved  to  be  very  useful  in  elucidating  pos¬ 
sible  mechanisms  of  Fe  embrittlement One  of  the  major  obstacles  in  using  the  first  principles 
approach,  however,  has  been  the  need  for  substantial  computing  power  not  available  until  recently. 
With  the  rapid  development  of  computer  technology,  one  can  now  perform  such  calculations  with 
reasonable  cost. 

In  a  previous  paper,  we  applied  the  first  principles  DMol  molecular  cluster  method  to  study 
phosphorus  induced  structural  changes  in  the  Fe  grain  boundary  (GB).'°  It  was  shown  that  P 
induces  large  relaxations  in  the  GB,  and  the  electronic  structure  of  the  GB  also  changes  significantly. 
In  this  paper,  we  compare  the  effects  of  B  and  S  impurities  on  the  Fe  S3[li0](lll)  grain  boundary 
using  the  same  theoretical  techniques:  We  first  determine  the  atomic  relaxation  induced  by  the 
impurities  using  atomic  force  calculations;  the  electronic  structure  is  then  calculated  for  the  relaxed 
structure.  Previous  electronic  structure  calculations  used  relatively  small  clusters^  and  did  not 
consider  the  role  of  lattice  relaxation.®’®  Those  studies  provided  useful  information  on  the  bonding 
between  the  metal  and  impurities.  However,  as  the  large  relaxations  induced  by  impurities  may 
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also  change  the  electronic  structure  of  the  Fe  GB,  consideration  of  atomic  relaxation  is  necessary 
for  a  correct  description  of  the  electronic  structure. 

II.  Methodology 

We  use  the  DMol  method^^  which  is  a  first-principles  numerical  method  for  solving  the  local 
density  functional  equations  and  is  capable  of  calculating  analytic  energy  gradients^^  for  each  atom 
within  the  cluster  model.  The  Hedin-Lundqvist  exchange  correlation  potential^^  is  employed  with 
the  frozen-core  approximation.  We  use  extended  basis  sets  for  the  B,  S  and  Fe  atoms,  ie.,  a  double 
set  of  valence  functions  plus  a  single  d  polarization  function  for  B  and  S  atoms  and  a  single  p 
function  for  Fe.  The  binding  energy  of  a  cluster  is  defined  as  Eb  =  -{Et  -  Ea)  where  Et  is  the 
total  energy  of  the  cluster  and  Ea  is  the  sum  of  atomic  total  energies.  For  a  given  atomic  geometry, 
the  binding  energy  of  the  system  and  the  forces  on  atoms  of  interest  are  calculated.  To  find  the 
optimized  geometry,  the  atoms  are  further  displaced  according  to  the  forces  acting  on  them.  An 
optimized  structure  is  obtained  when  all  the  forces  acting  on  the  atoms  are  sufficiently  small.  In 
this  work,  the  degree  of  self-consistent  convergence  is  measured  by  root  mean  square  (rms)  changes 
in  the  charge  density;  it  is  set  to  10"®  which  allows  the  total  energy  to  converge  to  10”=’  Ry.  The 
force  convergence  criterion  is  6.0  X  10"^  Ry/a.u.. 

The  grain  boundary  considered  is  the  E3[ll0](lll)  ’’incoherent”  twin  boundary  as  treated  in 
previous  studies.*"^®’^”*  The  91  atom  cluster  model  contains  17  layers  of  Fe  adjacent  to  the  grain 
boundary  and  is  shown  in  Fig.l.  The  impurity  atoms  are  placed  at  the  center  of  the  trigonal  prism 
formed  by  iron  atoms  in  the  GB  core,  and  D3/1  symmetry  is  used.  In  searching  for  the  optimum 
geometry,  the  boundary  atoms  of  the  cluster  are  fixed  in  their  bulk  positions.  When  we  refer  to 
the  relaxation  of  a  particular  atom,  say,  Fen,  other  atoms  satisfying  D3/1  group  operations  have 
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the  same  relaxation  as  Fen- 


III.  Results  and  Discussion 
A.  Atomic  relaxation  study 

The  relaxation  study  is  carried  out  by  force  calculations  in  which  atoms  in  columns  1  and 
2  in  Fig.  1  can  relax  along  the  2  direction  since  those  atoms  are  most  affected  by  the  impurity 
in  the  GB  center.  For  the  B/Fe  system,  we  first  calculate  the  relaxation  of  Fe  atoms  using  a  53 
atom  cluster  which  takes  less  computing  time.  Since  the  interlayer  distance  and  the  bond  length 
between  impurity  and  Fe  atoms  are  underestimated  using  the  53  atoms  cluster  as  previously  shown 
in  the  P/Fe  system, the  result  from  the  small  cluster  is  used  as  an  initial  input  for  the  91  atoms 
cluster  calculation.  The  final  optimized  structural  information  is  presented  in  Table  I.  For  the  S/Fe 
system,  we  use  the  P/Fe  system  data  for  the  91  atom  cluster^®  as  initial  input,  anticipating  that  S 
and  P  induce  similar  relaxations. 

The  basic  relaxation  results  are:  1)  for  the  clean  GB  structure,  there  are  two  structures  having 
energy  minima.  One  corresponds  to  the  case  in  which  two  center  Fe  atoms  (Fe2i)  move  towards 
the  GB  center,  referred  to  as  GB-1.  This  gives  a  structure  for  the  E3  GB  core  which  is  very  similar 
to  the  w  phase  observed  in  Ti  alloys.  The  other  structure  corresponds  to  Fe2i  atoms  moving  apart 
from  the  GB  center  (an  ’’anti-w”  structure),  referred  to  as  GB-2.  It  has  been  shown  that  the  GB-1 
structure  is  more  stable  energeticaUy  than  the  GB-2  structure.^^  Here,  we  also  allow  the  Fe32  atoms 
to  relax  along  the  .r  direction  in  the  GB-1  structure.  It  is  found  that  the  calculated  results  show 
an  average  difference  of  less  than  1%  for  structural  data  and  very  small  changes  for  binding  energy 
compared  to  the  GB-1  structure  displayed  in  Table  I.  Hence,  the  extra  relaxation  of  the  Fe32  atoms 
causes  oidy  slight  changes  for  the  clean  GB  structure.  2)  When  impurities  are  introduced  in  the 
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center  of  the  GB,  the  Fe  atoms  directly  above  the  impurity  (Fe2i  atoms)  are  pushed  outwards 
from  the  center  (favoring  the  GB-2  structure),  and  this  affects  the  atoms  in  the  outer  layer.  The 
net  effect  is  that  almost  all  Fe  atoms  are  relaxed  outwards  from  the  center.  3)  S  induces  larger 
relaxation  for  atoms  near  the  impurity  (eg.,  dn  and  d^i)  than  does  B  as  shown  in  Table  I.  The  Fe 
atom  (Feji)  directly  above  the  impurity  relaxes  outwards  42%  for  S/Fe  and  30%  for  B/Fe  compared 
with  its  initial  (unrelaxed  coincidence  site  lattice)  position.  The  bond  length  between  Fe2i  and  B 
is  2.15A  and  is  smaller  than  that  between  Fe2i  and  S  (2.35A). 

The  relcLxation  results  can  be  compared  with  an  empirical  calculation  using  a  modified  Em¬ 
bedded  Atom  Method  (EAM)  by  Krasko. The  EAM  calculations  show  that  for  the  clean  GB,  the 
stable  structure  corresponds  to  the  down  position  of  Fe2i  which  agrees  with  our  GB-1  structure. 
After  B  and  S  are  put  in  the  center  of  the  GB,  the  EAM  also  found  that  S  induces  larger  relax¬ 
ation  of  Fe  atoms  near  the  impurity  than  does  B.  The  calculated  bond  length  between  Fe2i  and  B 
and  between  Fe2i  and  S  from  the  EAM  are  1.96A  and  2.60A  respectively,  which  has  an  average 
difference  of  10%  with  the  results  of  the  present  calculation.  Our  bond  length  value  for  Fe2i  and  S 
(2.35A)  may  increase  if  an  even  larger  cluster  is  used.  As  can  be  seen  from  Table  I,  the  distance 
between  the  boundary  atoms  and  the  nearest  inner  atoms  in  the  0  direction  are  compressed  (ie.  di3 
and  d23)  compared  to  the  bulk  value  of  2.48A.  It  is  quite  possible  that  the  boundary  atoms  wiU  be 
pushed  further  outwards  when  a  larger  cluster  is  used.  If  that  is  true,  we  can  expect  a  weakening 
of  the  vertical  interaction  between  Fe2i  and  S.  This  may  be  useful  in  elucidating  the  mechanism  of 
S  embrittlement  since  the  competition  between  bonding  in  the  vertical  and  horizontal  directions  is 
an  important  factor. 

It  is  worth  mentioning  that  although  B  is  a  light  element,  it  has  a  higher  binding  energy 
(1.71  eV)  in  the  Fe  GB  than  does  S.  This  is  due  to  the  shorter  Fe2i  and  B  bond  length  that  results 
in  the  stronger  interaction  between  B  and  Fe.  Similar  results  were  obtained  for  the  binding  energy 
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of  B  and  S  in  Ni  by  Painter  and  Averill®  which  showed  that  B  binds  in  the  Ni  host  much  more 
strongly  (7.0  eV)  than  does  S  (2.2  eV)  and  produces  a  smaller  expansion  of  the  host  cluster  than 

does  S. 

B.  Charge  density  and  density  of  states 

Figure  2  shows  the  charge  density  difference  plotted  in  the  (lIO)  plane  containing  the  impurity 
and  several  neighboring  Fe  atoms.  (The  charge  density  difference  is  calculated  by  subtracting  from 
the  charge  density  of  the  impurity  in  the  GB,  the  pure  GB  charge  density  and  a  single  impurity 
charge  density  placed  at  the  center  of  the  GB  so  that  the  impurity-induced  charge  redistribution  can 
be  seen  more  clearly.)  Comparing  the  charge  density  difference  between  the  B/Fe  and  S/Fe  systems, 
we  can  see  that  there  is  one  common  feature  for  the  two  systems;  the  interaction  of  impurity  and 
Fe  atoms  is  restricted  to  a  local  region  near  the  impurity,  although  the  lattice  relaxation  is  known 
to  extend  beyond  this  area.^”* 

The  main  difference  between  them  is  the  charge  distribution  in  the  vertical  (2)  and  horizontal 
(i)  directions.  In  Fig. 2(a),  B  is  seen  to  form  relatively  strong  covalent-like  bonding  states  with  Fe 
atoms  in  the  vertical  direction  and  a  much  weaker  bonding  with  Fe  atoms  in  the  horizontal  direction. 
In  the  S/Fe  system  (see  Fig.2(b)),  sUghtly  more  charge  is  accumulated  in  the  region  between  S 
and  Feai  atoms  than  between  S  and  Fe2i  atoms,  giving  the  impression  that  the  bonding  along 
the  horizontal  direction  is  stronger  than  that  along  the  vertical  direction,  although  the  distance 
between  S  and  the  Fe  atoms  in  both  directions  are  almost  equal.  The  same  phenomenon  was 
found  for  the  P /Fe  system  calculated  both  from  the  full-potential  linearized  augmented  plane  wave 
(FLAPW^)  and  DMoP°  method.  It  is  doubtful  that  these  charge  accumulations  represent  strong 
covalent  bonding  between  S  and  Fe  since  the  density  of  states  (DOS)  curves  presented  below  do 
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FIG.  2.  Charge  density  difference  between  impurity  in  the  GB  and  a  superposition  of  a  single 
impurity  and  pure  GB  calculated  from  the  91  atom  cluster.  The  plots  are  in  a  plane  contaimng 
the  impurity,  Fen,  Feji  and  Feai  atoms.  The  contour  spacings  are  0.002e/(a.u.f .  Solid  Unes 
mean  again  of  charge,  and  dashed  lines  mean  a  loss  of  charge,  (a)  B/Fe. 


not  confirm  the  hybridization  between  S  and  Fe. 

To  understand  this,  we  plot  the  two  molecular  orbitals  containing  mainly  S-3p  electron  con¬ 
tributions  in  Fig.3(a)  and  Fig.3(b).  In  Fig.3(a),  one  can  see  a  dominant  S-3p^-Uke  wave  function 
having  only  weak  bonding  with  Feai;  a  similar  dominant  S-3px-like  wave  function  is  seen  in  Fig.3(b). 
Hence,  the  bonding  between  S  and  Fe  is  weak.  If  we  compare  Fig. 3  with  the  free  atom  S-3p  wave 
function  (not  shown  here),  we  find  that  the  free  S-3p  wave  function  has  greater  spatial  expansion. 
When  S  is  in  the  Fe  GB  environment,  the  S-3p  electrons  will  feel  the  repulsive  interaction  of  nearby 
Fe  electrons  and  shrink  towards  the  S  nucleus,  and  so  there  is  some  charge  accumulation  in  the 
region  between  the  Fe  and  S  atoms. 

Note  that  the  charge  accumulation  in  the  region  between  S  and  Fe  is  acompanied  by  a  decrease 
of  charge  in  the  S  center  region  as  can  be  seen  from  Fig.2(b).  To  estimate  how  the  charges  are 
distributed,  we  calculate  the '^Mulliken  electron  population;  the  results  show  that  there  is  about 
0.37e  charge  redistribution  from  S-3s  and  S-3p  states  to  the  high  energy  S-3d  state.  This  wiU 
reduce  the  charge  in  the  center  region  compared  to  the  free  S  atom  which  results  in  the  image  of 
Fig.2(b). 

Figures  4  and  5  show  the  calculated  partial  density  of  states  (DOS)  for  B,  S  and  several  nearby 
Fe  atoms.  The  DOS  of  the  clean  GB  (i.e.,  GB-1)  is  shown-  by  dashed  lines  and  the  DOS  of  single 
B  and  S  are  also  shown  by  dashed  lines  in  Figs.  4a  and  5a,  respectively.  For  S  in  the  Fe  GB,  the 
main  S-3s  peak  (not  shown)  lies  well  below  the  d  states  of  Fe,  and  thus  has  no  hybridization  with 
Fe.  The  majority  of  the  S-3p  states  and  a  small  portion  of  S-3s  states  are  at  the  tail  end  of  the 
Fe  states,  again  showing  a  weak  interaction  with  Fe.  For  B  in  the  Fe  GB,  the  B-2s  electron  has 
v0xy  little  overlap  with  Fe,  while  B"2p  forms  broadly  hybridized  states  with  Fe  d  states  between 
-2  to  -5.5  eV.  This  is  in  agreement  with  the  charge  density  picture  shown  in  Fig.2(a).  To  obtain 
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FIG.  4.  The  partial  density  of  states  for  S  in  the  Fe  GB. 
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fig.  5.  The  partial  density  of  states  for  B  in  the  Fe  GB. 
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FIG.  6. 


Contour  plot  of  the  loW  molecular  orbital  tor  B  in  the  Fe  GB.  The  contour  spacing  is 


0.004 


insight  into  how  B  bonds  with  Fe,  we  plot  in  Fig.  6  the  10a2"  orbital  which  is  located  2.82  eV 
below  Ef-  It  is  clearly  seen  that  B-2pz  and  Fe-3d^2  form  a  a  bonding  state  in  the  vertical  direction 
in  contrast  with  Fig.3(a)  which  shows  that  the  S-3p.  states  dominate  the  area  between  S  and  Fe 
without  overlapping  with  the  Fe  states. 

In  summary,  we  have  performed  first  principles  local  density  DMol  molecular  cluster  calcula¬ 
tions  on  the  role  of  the  B  and  S  impurities  in  the  Fe  grain  boundary.  When  impurities  are  placed  in 
the  center  of  the  GB,  it  is  found  that  B  induces  less  relaxation  for  nearby  Fe  atoms  than  does  S.  The 
nearest  distance  between  B  and  Fe  and  between  S  and  Fe  are  2.15A  and  2.35A,  respectively.  For 
the  bonding  between  impurities  and  Fe,  we  found  that  B-2p,  forms  a  relatively  strong  cr  bonding 
with  Fe2i-3pj2,  while  S  has  very  little  hybridization  with  Fe,  resulting  in  a  weak  bonding  state. 
These  differences  in  electronic  structure  for  B  and  S  may  be  useful  in  understanding  their  different 
effects  on  grain  boundary  cohesion. 

Acknowledgments 

This  work  was  supported  by  the  Office  of  Naval  Research  (Grant  No.  N00014-90-J-1363),  and 
by  a  computing  grant  at  the  Pittsburgh  Supercomputing  Center,  supported  by  the  NSF  Division 


for  Advanced  Scientific  Computing. 


References 


iG.B.  Olson,  in  Innovations  in  Ultrahigh-strength  Steel  Technology,  ed.  G.B.  Olson,  M.  Azrin  and 
E.S.  Wright,  Sagamore  Army  Material  Research  Conference  Proceedings;  34th(1990),  and 

references  therein. 

^J.R.  Rice  and  J-S.  Wang,  Mat.  Sci.&:  Eng.,  A  107,  23(1989). 

^R.P.  Messmer,  Phy.Rev.  B  23,  1616(1981). 

j^0ssmer  and  C.L.  Briant,  Acta  Metall.  30,  457(1982). 

^A.  CoUins,  R.C.  O’Handley  and  K.H.  Johnson,  Phys.  Rev.  B  38,  3665(1988). 

®G.S.  Painter  and  F.W.  Averill,  Phy.  Rev.  Lett.  58,  234(1987). 

’’M.E.  Eberhart  and  D.D.  Vvedensky,  Phys.  Rev.  lett.  58,  61(1987). 

®G.L.  Krasko  and  G.B.  Olson,  Solid  State  Comm.  76,  247(1990). 

®R.  Wu,  A.J.  Freeman  and  G.B.  Olson,  J.  Mater.  Res. 7,  2403  (1992). 

^•^S.  Tang,  A.J.  Freeman  and  G.B.  Olson,  Phys.  Rev.  B  47,  2441(1993). 
i^B.  DeUey,  J.  Chem.  Phys.  92,  508  (1990). 

12b.  DeUey  J.  Chem.  Phys.  94,7245(1991). 

i^L.  Hedin  and  B.I.  Lundqvist,  J.  Phys.  C4,  2064  (1971). 

i‘‘G.L.  Krasko,  in  Structure  and  Properties  of  Interfaces  in  Materials,  edited  by  W.A.T.  Clark,  U. 
Dahmen  and  C.L.  Briant,  MRS  Symposia  Proceedings  No.  238  (Materials  Research  Society, 
Pittsburgh,  1991),  p.  481 


10 


Tables 


TABLE  I.  Calcvilated  structural  parameters  (in  A)  for  the  91  atom  cluster  model  defined  in  Fig.l. 
The  ’’initial”  distances  are  calculated  from  the  unrelated  coincidence  site  lattice  structure.  GB-1 
refers  to  the  geometry  when  Fe2i  relaxes  inward  towards  the  center  of  the  GB  and  GB-B  and 
GB-S  refer  to  the  cases  when  B  and  S  are  added  to  the  GB;  dpt  is  the  distance  between  the  two 
near-center  atoms  in  column  2  in  the  pure  GB  and  Et,  is  the  binding  energy  in  eV. 


initial 

GB-1 

GB-B 

GB-S 

dll 

1.65 

2.49 

2.48 

2.61 

d\2 

2.48 

2.47 

2.41 

2.35 

d\z 

2.48 

2.07 

2.14 

2.14 

^21 

— 

— 

2.15 

2.35 

d22 

2.48 

2.48 

2.46 

2.30 

d23 

2.48 

2.93 

2.01 

1.97 

dpe 

3.31 

2.43 

— 

— 

Ek 

— 

421.86 

428.67 

426.96 
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